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Abstract. We present an automated tool for measuring atmospheric parameters (Teff, 
log g), [Fe/H]) for F-G-K dwarf and giant stars. The tool, called GAUFRE, is written 
in C++ and composed of several routines: GAUFRE-RV measures radial velocity from 
spectra via cross-correlation against a synthetic template, GAUFRE-EW measures atmo- 
spheric parameters through the classic line-by-line technique and GAUFRE-CHI2 per- 
forms ax^ fitting to a library of synthetic spectra. A set of F-G-K stars extensively studied 
in the literature were used as a benchmark for the program: their high signal-to-noise and 
high resolution spectra were analyzed by using GAUFRE and results were compared with 
those present in literature. The tool is also implemented in order to perform the spectral 
analysis after fixing the surface gravity (log g) to the accurate value provided by astero- 
seismology. A set of CoRoT stars, belonging to LRcOl and LRaOl fields was used for first 
testing the performances and the behavior of the program when using the seismic log g. 

1 Introduction 

Spectroscopy is one of the most powerful tool that astronomy possesses in order to derive atmospheric 
parameters and abundances of stars. From the stellar spectrum it is possible to measure the effective 
temperature (Teft), the surface gravity (log g) and the abundance of iron ([Fe/H]) and of several other 
elements. 

The classic method consists in measuring the equivalent widths (EW) of a species in two different 
ionization states, usually Fel and Fell. By imposing excitation and ionization equilibrium through 
stellar atmosphere models, it is possible to derive Tgff, log g and to infer elemental abundances from 
the curves of growth. This method is precise and it is widely adopted [1 J. In spite of the precision, the 
line-by-line classical analysis is time consuming: usually EWs are measured by hand using the IRAF 
"splot" routine (where the choice of the continuum and the position of the line is completely manual) 
and the procedure for finding the best Teft and log g requires a large number of iterations. 
The growing amount of stellar data due to large surveys (i.e. RAVE, SEGUE, LAMOST, HERMES 
and Gaia ESO) and the availability of dedicated telescopes and multi-object spectrograph, requires 
the development of automated pipelines, methods and programs that fasten the analysis process. For 
example, DAOSPEC [2| and ARES |3| are useful and free codes that performs an automated EW 
measurement; the MOOG code [41 is a valid collection of routines useful for determining atmospheric 
parameters and abundances. In particular, MOOG abfind and synth tasks are widely used in literature 

an]. 

Another method for the estimation of atmospheric parameters and elemental abundances is to compute 



e-mail: valentini@astro.ulg.ac.be 



EPJ Web of Conferences 



Table 1. Comparison between the radial velocity we measured with GAUFRE-RV (adopting the library of syn- 
thetic spectra of Fossati) and lAU values. All radial velocities are in km s"' 



Star 


Instrument 


lAU 


GAUFRE-RV 










Vrarf 




HD 003712 


Asiago-Echelle 


-3.9 


0.1 


-4.1 


0.3 


HD 012929 


Asiago-Echelle 


-14.3 


0.2 


-13.8 


0.4 




ESO-FEROS 






-14.6 


0.2 


HD 062509 


Asiago-Echelle 


-1-3.3 


0.1 


+3.3 


0.6 




ESO-UVES 






+3.4 


0.3 


HD 065934 


Asiago-Echelle 


+35.0 


0.3 


+34.6 


0.5 


HD 090861 


Asiago-Echelle 


+36.3 


0.4 


+36.7 


0.4 


HD 212943 


Asiago-Echelle 


+54.3 


0.3 


+54.2 


0.4 




ESO-UVES 






+54.1 


0.2 


HD 213014 


Asiago-Echelle 


-39.7 


0.0 


-40.0 


0.5 



a set of synthetic spectra and to find the best match between the synthetic and the observed spectrum. 
This method is adopted by surveys as RAVE |6| and HERMES [7| as well as in small surveys as ARCS 
|8| |9|. 

In this paper we present a new automatic code, GAUFRE, that can perform both type of analysis. In 
section[2]we give a short description of the idea behind the program and the two types of analysis that 
it performs. In subsection |2.4| we present the first results of tests using spectra of objects well known 
in literature. In section [3] we present an additional tool of GAUFRE that consists in using the astero- 
seismic gravity as a fixed value for log{g) in order to refine the measurement of Tefj , microturbulence 
velocity {^mic) ^nd elemental abundances. Section |4]discusses the future perspectives for the code. 



2 The GAUFRE code 



GAUFRE is a collection of several C++ routines. It is written in order to measure fast and precisely 
radial velocity (Vrad) and atmospheric parameters (Teff, log g, [Fe/H]) of a star starting from a one- 
dimensional normalized spectrum. The radial velocity is measured with a cross-coiTelation technique 
(routine GAUFRE-RV), atmospheric parameters can be measured adopting a fitting over a library 
of synthetic spectra (GAUFRE-CHI2 routine) or with the classic Fel-Fell lines technique (GAUFRE- 
EW routine). 

GAUFRE is an updated and extended version of the code written for the Asiago Red Clump Spectro- 
scopic survey (ARCS) |8|. GAUFRE was created because of the need of introducing new libraiies of 
synthetic spectra, adapting the code for different resolutions and implementing the classic line-by-line 
technique. The code was written in C++ and does not need any particular library, in order to avoid 
any software license problems and to be executable in different platforms. The only additional pro- 
gram needed for running GAUFRE is MOOG, which can be easily installed and downloaded from 
its homepage (http://www.as.utexas.edu/~chris/moog.html). The user is also supposed to download 
the libraries of synthetic spectra and model atmospheres required by GAUFRE. So far GAUFRE is 
tested for F-G-K stars (both giants and dwarfs) and works in the 3500 - 9000 A spectral range. 
The program has currently been tested for spectral resolutions of the Asiago Echelle Spectrograph 
(R=20,000), ESO-FLAMES-GIRAFFE HRIO, HR15N, HR21, HR14 setups (R=16,000 - 20,000), 
ESO-FLAMES-UVES U580 setup (R=52,000), ESO-FEROS (R=48,000) and AAOmega 2dF 580V 



and 385R (R= 1,300). Details about the adopted libraries can be found on subsections 2.3 and 2.4 
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Table 2. Libraries of synthetic spectra used by GAUFRE-RV and GAUFRE-CHI2. 



Library 


Model 


Linelist 


Solar abundances 


Resolution 




atmospheres 








Munari et al. (2005) 


Castelli &. Kurucz (2003) 


Kurucz 


Grevesse & Sauval (1998) 


20,000 


Fossati 


Castelli & Kurucz (2003) 


VALD2 


Grevesse & Sauval (1998) 


1,000,000 


De Lavemy et al. (2012) 


MARCS 




Grevesse, Asplund & Sauval (2007) 


150,000 



2.1 Radial Velocities: GAUFRE-RV 

The radial velocity subroutine (GAUFRE-RV) measures the radial velocity Viad by cross-correlating 
the observed spectrum with a synthetic spectral library [iO|. The procedure is the same as described 
in IS]. The procedure starts from the continuum normalized spectrum in a 2-column ASCII format. 
First, the synthetic normalized spectra of the selected library are renormalized, following the same 
parameters as adopted for the normalization of the observed spectrum (same function, same order and 
same high and low rejection values). To lower the impact of the different noise level in the observed 
and synthetic spectra, they are scaled to match their geometric mean. This procedure is needed in or- 
der to improve the accuracy of the cross-correlation and of the fitting and it is performed by the 
GAUFRE-LIB routine (see also section 2.2). 

The result of the GAUFRE-RV subroutine is a file containing the value of Vjad and an ascii file con- 
taining the continuum normalized spectrum corrected for the radial velocity. 

To validate the GAUFRE-RV routine we measured the radial velocity of a set of Red Giants stars that 
are lAU standard radial velocity stars. These stars were observed with the Asiago Echelle Spectro- 
graph (INAF-OAPd); we then downloaded, when available, the spectrum from the ESO- Archive. For 



this test, we adopted the synthetic spectra library provided by L. Fossati (cf. section 2.2 1. Results are 
summarized in Table [T] The mean difference between our values and those present in the literature is 
ARVq = 0.10 km s^', with arms of 0.3 km s^'. 



2.2 x"^ on synthetic spectra libraries: GAUFRE-CHI2 

Atmospheric parameters (Teff , log g, [Fe/H], [a/Fe], Vjot sin /) are obtained by GAUFRE-CHI2 wiax^ 
fitting of the continuum normalized spectrum against a synthetic spectral library. 
The choice of the spectral library is up to the user, available libraries are: a library based on Kurucz 
model atmosperes |11 1, the library provided by L. Fossati (built adopting the spectral synthesis code 
Synth3 described in 1 12|) and the AMBRE library 1 13 1. Characteristics of different libraries are sum- 
marized in Table [3] Libraries are provided at different resolutions and they cover a wide spectral range 
(usually 3,500 - 10,000 A). Before starting x^ analysis, the desired library must be cut, normalized 
and degraded at the same wavelength interval, normalization function and resolution of the real spec- 
tra. For this purpose a routine GAUFRE-LIB has been created. The degradation of the spectra at the 
desired resolution is performed through deconvolution with a Gaussian profile. 



2.3 EWand MOOG: GAUFRE-EW 

GAUFRE-EW automatically performs the classical line-by-line analysis for deriving atmospheric pa- 
rameters and abundances. The procedure starts from a continuum normalized spectrum (ASCII 2 
columns format), a file containing a list of lines to measure and their parameters (as requested by 
MOOG) and a file containing the parameters of the spectrum like the wavelength coverage, resolution 
and, if any, guessed values for Ten and log g. 

The EW of every line present in the input file is measured (except when it is not detectable). The 
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Fig. 1. Left panel: Comparison between the EW measured by GAUFRE-EW (y-axis) and the EW measured with 
the IRAF-splot routine, fitting the line to a Gaussian profile. Right panel: Differences between the values of 
(top panel), log g (middle panel) and [Fe/H] (bottom panel) measured by GAUFE-EW and GAUFRE-CHI2 and 
those present in literature. 



program selects an area of 3-4 A around the wavelength of the line (this parameter is selected by 
the user). The spectrum is then fitted with a polynomial function in order to determine the continuum 
and point with the lowest intensity. In order to test the automated EW measurement of the line we 
compared the EW values obtained by measuring features by hand (using IRAF-splot) and the corre- 
spondent EW measured by GAUFRE. The test was performed on the spectrum of Arcturus, taken with 
the ESO-FLAMES-UVES instrument^ith the U580 setup (5770-6825 A). The agreement is quite 
good and rms is of ~3 mA (see Figurellb. 

Atmospheric parameters were computed by using MOOG abfind driver (the program uses its non- 
iinteractive version, MOOGSILENT), the measured EW of Fel and Fell lines and a family of model 
atmospheres (MARCS 1 14| or Kurucz 1 15|). Teff is calculated by assuming the excitation equilibrium 
and minimizing the trend of the Fe abundance versus the excitation potential. The surface gravity, log g 
is derived by assuming the ionization equilibrium: log n (Fel)= log « (Fell). The procedure is iterative 
and the program will converge to Ten and log g that satisfy both the ionization and excitation equi- 
libria. The value of the microturbulence ^mic is derived by minimizing the trend of the Fel abundance 
versus the Fel lines EW. 



2.4 Atmospheric parameters validation 

We took spectra of 7 F-G-K stars taken with FLAMES-UVES U580 or FEROS from the ESO-archive 
(http://archive.eso.org ). We selected spectra of targets very well known in literature: a Cen A, jj Cas A, 
P Vir, Arcturus, /i Leo, ^ Hya and y Sge. As a reference we used an average value of the most recent 
entries of the PASTEL catalog 13 |. We analyzed spectra with both GAUFRE-EW (Kurucz model 
atmospheres) and GAUFRE -CHI2 (Fossati library of synthetic spectra) and we compared our results 
with those present in literature. The agreement is quite good as showed in Tab.3. 
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Table 3. Comparison of the atmospheric parameters obtained by GAUFRE-CHI2 (G-CHI2) and GAUFRE-EW 
(G-EW) with thopse present in literature for a set of 7 stars. 



GAUFRE-CHI2 






GAUFRE-EW 






(T'cfflo-CHK - Teff) = 16 K 

(log gr|G-cHi2 - log 9) = 0.13 dex 
<[Fe/H]|o-cHi2 - [Fe/H])=-0.05 dex 


cr = 
cr = 
cr = 


84 K 
0.17 dex 
0.10 dex 


(TefflG-EW - ^eff) = K 

(log gb-EW - log 9> = 0. 1 1 dex 
<[Fe/H]|G-EW - [Fe/H]> =-0.06 dex 


<T = 

a = 
(T = 


17 K 
0.17 dex 
0.10 dex 



3 Asteroseismic constraints 

As extensively discussed, for example, in lfT6l . the frequency of maximum power, Vmax can be used for 
deriving a precise estimate of the surface gravity: 

log g = log go + log I I + :^ log i:^^] ■ (1) 



2 \ Teff , / 

The precision of the gravity derived from Eq.[T]is generally expected to be below 0.05 dex, thanks to 
the weak sensitivity of the scaling relation to the assumed Tgff and the high precision usually achieved 
for the measurement of Vmax- 

The scaling relation has been demonstrated to be reliable ITTI ifTSl ^19\ pOl and it can be used for re- 
fining atmospheric parameters and abundances. As a matter of fact one can improve the spectroscopic 
analysis by fixing the log g to the seismic value: this largely increases the accuracy of the derived T^ff, 
[Fe/H] and hence chemical abundances. 

A set of spectra of 1 1 1 RG stars belonging to the LRcOl and LRaOl fields of CoRoT have been used as 
benchmark for testing the log g values derived by GAUFRE. Spectra are taken with the ESO-FLAMES 
GIRAFFE 9 setup, centered in the Mgl triplet (5278 A). Spectra have already been analyzed [211 by 
using MATISSE tool 1221 . 

First, we compared log g values of GAUFRE-CHI2 (using the synthetic library provided by L. Fos- 
sati, see Tat[3| with those provided by asteroseismology (see Fi^ top panels), then we compared the 
values of Teff and [Fe/H] derived by fixing log g to the seismic value (see Fi^ middle and bottom 
panels). The set of spectra taken from the work of Gazzano has very poor SNR (see Fi^2]i, of ~30 on 
average. In addiction the spectral range covered by spectra is affected by the presence of MgH molec- 
ular bands. This complicates even more the continuum normalization, leading to several systematics 
in the atmospheric parameters determination. 

As shown in Fig. [2j the seismic gravity can be used to enhance the accuracy of the atmospheric 
parameters. In the case of Gazzano et al. dataset |21 1 with a poor SNR, the seismic gravity helped in 
providing a more precise determination of Teff and [M/H] reducing average eiTors from 120 K and 
0.20 dex to 75 K and 0.11 dex respectively. 



4 Conclusions 

We present the GAUFRE program, a versatile tool developed for measuring radial velocities and at- 
mospheric parameters from optical spectra. 

The program is composed by different routines: GAUFRE-RV, GAUFRE-CHI2, GAUFRE-EW. 
We performed some preliminary tests in order to check the performances of the tool. These tests show 
that the program is reliable and that it can be used to process spectra of F-G-K dwarfs and giants with 
a SNR above 40. At the moment it is adopted by the Liege node within the Gaia-ESO Survey (PI: G. 
Gilmore and S. Randich) and further applications are planned. 

We also showed an interesting and useful extension of GAUFRE that uses, when avaiable, the seismic 
log g as a fixed value for the surface gravity: the precise and also likely more accurate values given by 
asteroseismology allow us to greatly refine the values of Tgff and [Fe/H]. 

The GAUFRE tool is continously in development: we plan to implement new libraries of synthetic 
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Fig. 2. Left panel: differences between the atmospheric parameters measured with different techniques (GAUFRE- 
CHI2 this paper, MATISSE | 21 1, photometry (J - K) and asteroseismology) in function of the SNR. Right panel: 
differences between the atmospheric parameters measured with different techniques (GAUFRE-CH12 this paper, 
MATISSE by |21 1, photometry (J — K) and asteroseismology) in function of the Teff, log g and [M/H] measured 
by GAUFRE-CHI2 by fixing the log g to the seismic value. Green crosses are data from |2l |, blue diamonds 
are data obtained by using GAUFRE-CHI2 and red circles are values measured by GAUFRE-CHI2 by fixing the 
log g to the seismic value. It is worth to take into account the poor quality of the spectra: more than 50% of the 
spectra possess a SNR < 40. 



spectra and to extend the analysis to the infrared region. New and detailed tests are planned as well, in 
order to better investigate the performances of GAUFRE at different SNR. 

In the next future we plan to make the GAUFRE code avaiable through the web and to create an user- 
friendly graphical interface. 
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